Agnafors et al. Child and Adolescent Psychiatry and Mental Health 201 3, 7:1 0 
http://www.capmh.eom/content/7/1/10 



CHILD & ADOLESCENT 
PSYCHIATRY & MENTAL HEALTH 



RESEARCH Open Access 



Effect of gene, environment and maternal 
depressive symptoms on pre-adolescence 
behavior problems - a longitudinal study 

Sara Agnafors 1 *, Erika Comasco 2 , Marie Bladh 3 , Gunilla Sydsjo 3 , Linda DeKeyser 3 , Lars Oreland 2 
and Carl Goran Svedin 1 



Abstract 

Background: Depression is a common and disabling condition with a high relapse frequency. Maternal mental 
health problems and experience of traumatic life events are known to increase the risk of behavior problems in 
children. Recently, genetic factors, in particular gene-by-environment interaction models, have been implicated to 
explain depressive etiology. However, results are inconclusive. 

Methods: Study participants were members of the SESBiC-study. A total of 889 mothers and their children were 
followed during the child's age of 3 months to 12 years. Information on maternal depressive symptoms was 
gathered postpartum and at a 12 year follow-up. Mothers reported on child behavior and traumatic life events 
experienced by the child at age 12. Saliva samples were obtained from children for analysis of 5-HTTLPR and BDNF 
Val66Met polymorphisms. 

Results: Multivariate analysis showed a significant association between maternal symptoms of depression and 
anxiety, and internalizing problems in 12-year-old children (OR 5.72, 95% CI 3.30-9.91). Furthermore, carriers of two 
short alleles (s/s) of the 5-HTTLPR showed a more than 4-fold increased risk of internalizing problems at age 12 
compared to l/l carriers (OR 4.73, 95% CI 2.14-10.48). No gene-by-environment interaction was found and neither 
depressive symptoms postpartum or traumatic experiences during childhood stayed significant in the final model. 

Conclusions: Concurrent maternal symptoms of depression and anxiety are significant risk factors for behavior 
problems in children, which need to be taken into account in clinical practice. Furthermore, we found a main effect of 
5-HTTLPR on internalizing symptoms in 12-year-old children, a finding that needs to be confirmed in future studies. 
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Background 

Depression and anxiety conditions are an increasing 
problem, leading to substantial economic and social con- 
sequences [1]. The prevalence of depression in pre- 
adolescents is approximately 1% [2] and there is an 
increase in rates of various psychiatric problems during 
adolescence and adulthood in individuals with a history 
of childhood depression [3,4]. Moreover, childhood de- 
pression has shown a risk factor pattern different from 
that of adolescent and adult onset depression, raising the 



* Correspondence: sara.agnafors@liu.se 

1 Division of Child and Adolescent Psychiatry, IKE, Faculty of Health Sciences, 

Linkoping University, S-581 85 Linkoping, Sweden 

Full list of author information is available at the end of the article 



question whether childhood depression is etiologically 
separate from the latter [5,6]. 

The longitudinal approach serves as an excellent way of 
studying risk factors, development and course of psychi- 
atric disorders and symptoms. Many factors may elevate 
the risk for onset of depression and emotional problems 
in children. Postpartum depression has been shown to 
predict child behavior problems in numerous studies [7,8], 
but there is also support for the fact that on-going mater- 
nal depression exerts a risk factor for internalizing as well 
as externalizing symptoms in children [9-11]. 

Among other known risk factors for depression and emo- 
tional problems are traumatic life events [12] - children ex- 
posed to physical abuse, parental divorce and domestic 
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violence are at increased risk of behavior problems and 
adult depression [13-15]. Still, not everyone subjected to 
traumatic life events develops depressive symptoms. Hence, 
there may be individuals predisposed or more vulnerable to 
certain risk factors. 

In recent years, much effort has been given to explaining 
depression in relation to genetic factors. Many candidate 
genes have been examined, of which the serotonin trans- 
porter gene (S-HTT) is the most studied. The 5-HTT gene 
includes a functional polymorphism, the serotonin trans- 
porter gene-linked polymorphic region (5-HTTLPR), 
which consists of two common alleles; short (s) and long 
(1). Less effective serotonin expression and availability have 
been shown in s-allele carriers, compared to individuals 
homozygous for the 1 allele [16]. The gene-by- environ- 
ment model put up by Caspi et al [17], showing an associ- 
ation between the 5-HTTLPR, stressful life events and 
depression, has since publication been subjected to nu- 
merous attempts at replication. Some researchers have 
been able to reproduce the results partially or completely 
[18-24], whereas others have not [25-28], however 
between-study heterogeneity has to be taken into account. 
When looking particularly at 5-HTTLPR gene-by-environ- 
ment studies on children, research is not that extensive 
and results have been conflicting. Araya et al [7] found no 
association between 5-HTTLPR and emotional symptoms 
in seven-year-olds whereas associations have been shown 
by other researchers, although on smaller study samples 
[20,22,29,30] . Moreover, also direct association between 5- 
HTTLPR and depression has been found [16,31-34]. 

Given not only the inconsistency of single genetic ex- 
planatory models, but also the complex etiology of depres- 
sion, researchers have lately been looking at gene-by-gene 
(-by-environment) models for an explanation of depressive 
vulnerability. Another candidate gene for depressive disor- 
ders is the Brain Derived Neurotrophic Factor (BDNF), 
which has previously been shown to act in synergy with 
serotonin [35]. BDNF is involved in reparation, plasticity 
and neurogenesis in the brain, and a single nucleotide poly- 
morphism (SNP) G/A (Val66Met) in the BDNF gene has 
been shown to affect levels of BDNF in the brain [36]. Simi- 
larly to 5-HTTLPR, studies regarding the association be- 
tween BDNF Val66Met and depression have been both 
confirmative [37,38] and negating [39]. 

A candidate gene-by-gene-by-environment interaction ef- 
fect of 5-HTTLPR-by- J BDA/FVal66Met-by-childhood adver- 
sity on depression has been shown [40] and to our 
knowledge four studies have attempted replication of this 
three-way interaction effect, however with contradictory re- 
sults [41-44]. The findings are inconsistent regarding both 
the presence of a three-way interaction effect and the risk 
genetic variants in the presence of childhood adversity. 
Two studies found evidence of a three-way interaction 
effect [42,44] whereas two other studies did not [41,43]. 



In view of this, the aim of the present paper was to study 
the gene-by-environment interaction on depressive symp- 
toms in 12-year-old children. We hypothesized to find a 
gene-by-environment, and possibly gene-by-gene-by-envir- 
onment interaction on depressive symptoms, including the 
genetic polymorphisms 5-HTTLPR and BDNF Val66Met, 
traumatic life events and maternal symptoms of depression. 

Methods 

Population and procedures 

This study is part of the South East Sweden Birth 
Cohort study (SESBiC-study) which started in 1995 with 
the purpose of early identification of psychosocially bur- 
dened families where children were at risk of dysfunc- 
tional development. Follow-ups have been carried out at 
ages 3, 5.5 and 12 and have been reported previously 
[11,45-48]. 

Baseline 

All mothers of children from a birth cohort born be- 
tween May 1st 1995 and December 31st 1996 in south- 
ern Sweden were asked to take part in the study, of 
whom 1723 mothers (88%) agreed to participate. The 
mean age of the mothers was 28.2 ± 4.6 years at child- 
birth. Ninety six percent of the mothers (n=1574) were 
cohabitating, 3.5% (n=57) were single parents, while 
0.5% (n=8) reported other family arrangements. Most 
mothers were born in Sweden (88.6%), (n=1482), but 
6.2% (n=103) were born in Europe (excluding Sweden), 
and 5.3% (n=88) outside Europe. Of the newborn chil- 
dren, 52.8% were boys and there were 27 twin pairs. The 
baseline study was carried out at Child Welfare Centers, 
(CWC) in connection with the routine 3-month check- 
up. Questionnaires were administered and a psychologist 
also interviewed the mothers. 

12-year follow-up 

Current home addresses for all 1 723 families were 
obtained from the Swedish Tax Offices. An information 
letter and a consent form were sent to parents (i.e. legal 
guardians). Parents who did not return the consent form 
within three weeks were contacted by phone. A separate, 
simplified information letter was enclosed for the child. 
Two children and four mothers were deceased, ten had 
moved out of the country and 24 were learning disabled 
and could therefore not participate. These subjects were 
excluded from the original 1723 in the baseline study, 
which left 1687 eligible participants, of whom 889 
(52.7%) accepted participation. The follow-up was car- 
ried out at school where research assistants met with the 
children in small groups. The children provided saliva 
samples and filled out a package of questionnaires separ- 
ately (as part of a larger study). A package of question- 
naires was sent to the mothers' home addresses. Families 
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who had moved out of the original catchment area were 
contacted by mail and phone in accordance with the 
regular routine. Those who agreed to participate re- 
ceived questionnaires and saliva sampling kits by mail, 
or if they preferred, were visited by a research assistant 
and the survey was carried out at the child's home. 

Genetic analyses 

The non-invasive and all-in-one Oragene® DNA Collec- 
tion Kit (DNA Genotek) was used for the collection, 
stabilization and transportation of saliva samples. DNA 
was isolated according to the laboratory protocol for 
manual purification of DNA. 

BDNF Val66Met A/G SNP (rs6265) and 5-HTTLPR 
genotyping were carried out according to previously 
published protocols [49]. The genotyping was performed 
blind to psychosocial data. To estimate the quality-rate 
of genotyping errors, a random repetition of -13% of 
the sample was carried out; the comparison indicated no 
inconsistencies. The genotypes were in Hardy- Weinberg 
equilibrium (5-HTTLPR (/ =2.73; p=0.10); females 
(^=0.44; p=0.51); males (/=2.67; p=0.10), BDNF 
Val66Met (^=0.08; p=0.77); females (/=1.20; p=0.27); 
males (^=1.98; p=0.16)). 

Questionnaires 
Baseline 

The Edinburgh Postnatal Depression Scale (EPDS) [50] is 
a widely used self-report questionnaire designed to screen 
for post-natal depression. EPDS refers to the 7 days pre- 
ceding completion of the form and was filled out by the 
mothers at baseline. 

Life Stress Score (LSS) is a 50-item semi- structured 
interview form, consisting of three main domains regard- 
ing the mother's social situation, medical information and 
psychological circumstances. The LSS has been used pre- 
viously in a Swedish population based study [51] and was 
filled out by a psychologist after interviewing the mothers 
at baseline. 

12-year follow-up 

The Child Behavior Check List (CBCL) [52] is a 113- 
item form assessing child behavior, focusing on subscales 
of internalizing and externalizing behavior respectively. 
We used the CBCL 4-18 years, which was filled out by 
the mothers at the child's age of 12. 

The Hopkins Symptom Checklist (HSCL-25), a short 
version of the HSCL-90 [53] was used to measure symp- 
toms of anxiety and depression during the most recent 
14 days. The HSCL-25 was filled out by the mothers at 
the 12-year follow-up. 

To assess potentially traumatic life events experienced 
by the child, the Swedish version of Life Incidence of 
Traumatic Events (LITE) was used [54,55]. The 16-item 



parent report form (LITE-P) was filled out by the 
mothers at the child's age of 12. 

Data analysis 

In accordance with previous studies, we used two cut offs 
for the EPDS in order to catch depressive symptoms of dif- 
ferent severity. The first (cut off 10) representing a screen- 
ing level [50] and the second (cut off 13) supposed to catch 
more severe clinical depressive symptomatology [56]. For 
the HSCL-25 total score, a mean item score of 1.75 was 
used as cut off, as has been used previously [57]. On the 
CBCL and LSS scales, the 90th percentile was set as a cut 
off. Results on the LITE form were dichotomized into 0-2 
events and >3 events, for which there is support in the 
literature [17]. Bivariate analyses between genetic markers 
(5-HTTLPR and BDNF Val66Met) and psychological scales 
(CBCL), as well as between scales (EPDS, HSCL-25, 
CBCL), were performed using the chi-square statistic. 
Multivariate analyses, with CBCL scales as dependent vari- 
ables and psychological scales, socio-demographic variables 
(LSS) and genetic markers as independent variables were 
also performed. Ethnic background (both parents born in 
Sweden, compared to one or both parents born abroad) 
and sex of the child were also controlled for. The multivari- 
ate analysis consisted of conditional stepwise logistic regres- 
sion considering full factorial models. However, since this 
procedure may choose models containing interactions 
without corresponding main effects, the models have been 
corrected for this and further evaluated and reduced to in- 
clude models with significant main effects and appropriate 
corresponding interactions. Results are presented with cor- 
responding Odds Ratios (OR) and 95% Confidence Inter- 
vals (CI). All statistical analyses were performed using IBM 
SPSS version 19 (IBM Corporation, Armonk, NY). 

Dropout rate analysis 

At the 12-year follow-up, the total dropout rate was 47.3% 
(n=798). There was a difference when comparing immi- 
grant status between participants and non-participants at 
the 12 year follow-up, where 54.6% (n=802) of mothers 
born in Sweden (n=1468) took part compared to 44.6% 
(n=45) of mothers born in Europe (n=101) and 34.9% 
(n=30) of mothers born outside of Europe (n=86) 
(x2= 15.79, p=0.00). Likewise, differences were found be- 
tween participants and non-participants at the follow-up 
when comparison for socio-demographic factors at base- 
line was made. Of mothers who scored above cut-off on 
the LSS total scale at baseline (n=141), 60.3% (n=85) took 
part in the follow-up compared to 70.7% (n=1093) of 
mothers with low total score at baseline (n=1546) 
(x2=6.65, p=0.01). No differences were found between par- 
ticipants and non-participants at the 12-year follow-up re- 
garding symptoms of postpartum depression. 
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Table 1 Odds ratios in predicting CBCL subscales at age 12 



Odds ratio 



95.0% CI for odds ratio 



p-value 



Internalizing symptoms > 90th percentile 

BDNF Val/Met compared to Val/Val 
BDNF Met/Met compared to Val/Val 
5HTTLPR s/l compared to l/l 
5HTTLPR s/s compared to l/l 
LITE 

EPDS screening cut off 10 

Clinical cut off 13 

HSCL-25 

LSS 

Sex 

Ethnicity 

Anxious depressed > 90th percentile 

BDNF Val/Met compared to Val/Val 
BDNF Met/Met compared to Val/Val 
5HTTLPR s/l compared to l/l 
5HTTLPR s/s compared to l/l 
LITE 

EPDS screening cut off 10 

Clinical cut off 13 

HSCL-25 

LSS 

Sex 

Ethnicity 

Withdrawn depressed > 90th percentile 

BDNF Val/Met compared to Val/Val 
BDNF Met/Met compared to Val/Val 
5HTTLPR s/l compared to l/l 
5HTTLPR s/s compared to l/l 
LITE 

EPDS screening cut off 10 

Clinical cut off 13 

HSCL-25 

LSS 

Sex 

Ethnicity 

Externalizing symptoms > 90th percentile 

BDNF Val/Met compared to Val/Val 
BDNF Met/Met compared to Val/Val 
5-HTTLPR s/l compared to l/l 
5-HTTLPR s/s compared to l/l 
LITE 

EPDS screening cut off 10 
Clinical cut off 13 



1.37 
0.82 
1.78 
4.43 
1.70 
1.90 
3.09 
5.92 
4.03 
0.63 
1.93 

0.99 
0.64 
1.10 
2.32 
1.72 
1.84 
2.72 
4.76 
3.53 
0.59 
1.19 

1.53 
0.75 
1.27 
2.48 
1.32 
2.04 
2.68 
4.54 
1.52 
1.28 
2.22 

1.28 
1.50 
1.60 
2.40 
1.62 
2.07 
2.52 



0.82-2.29 
0.19-3.52 
0.84-3.76 
2.09-9.38 
1 .04-2.80 
0.98-3.68 
1 .30-7.34 
3.59-9.77 
1 .95-8.35 
0.38-1.04 
1 .02-3.67 

0.59-1.64 
0.15-2.76 
0.59-2.06 
1 .22-4.42 
1 .07-2.72 
0.97-3.49 
1.15-6.43 
2.95-7.70 
1.71-7.28 
0.37-0.96 
0.59-2-39 

0.95-2.51 
0.17-3.26 
0.68-2.38 
1 .29-4.77 
0.83-2.10 
1.09-3.81 
1.13-6.34 
2.82-7.30 
0.62-3.69 
0.81-2.03 

1.22- 4.01 

0.78-2.10 
0.51-4.46 
0.86-3.01 

1 .23- 4.75 
1 .02-2.56 
1.13-3.80 
1.07-5.97 



0.24 
0.78 
0.13 
<0.01 
0.04 
0.06 
0.01 
<0.001 
<0.001 
0.07 
0.04 

0.96 
0.55 
0.76 
0.01 
0.03 
0.06 
0.02 
<0.001 
0.001 
0.03 
0.63 

0.80 
0.70 
0.46 
0.07 
0.24 
0.03 
0.02 
<0.001 
0.36 
0.30 
0.01 

0.33 
0.47 
0.14 
0.01 
0.04 
0.02 
0.03 
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Table 1 Odds ratios in predicting CBCL subscales at age 12 (Continued) 



HSCL-25 


5.20 


3.26-8.28 


<0.001 


LSS 


2.07 


0.93-4.62 


0.08 


Sex 


1.86 


1.17-2.95 


0.01 


Ethnicity 


1.89 


1 .04-3.46 


0.04 



Note: CI = Confidence Interval, CBCL = Child Behavior Checklist, BDNF = Brain Derived Neurotrophic Factor, 5-HTTLPR = serotonin transporter gene-linked 
polymorphic region, LITE = Life Incidence of Traumatic Events, EPDS = Edinburgh Postnatal Depression Scale, HSCL-25 = Hopkins Symptom Checklist 25, 
LSS = Life Stress Score. 

Logistic regression. Dependent variables: CBCL (internalizing symptoms, internalizing symptoms divided into 1) anxious depressed and 2) withdrawn depressed, and 
externalizing symptoms). Independent variables (0 is used as a reference level): LITE (0< 2, 1>3), EPDS (0<9, 1>10) and (0<12, 1>13), HSCL-25 (0=mean score <1.75, 
1= mean score >1.75): LSS (0<90th percentile, 1>90th percentile) sex (0=girls and 1=boys) and ethnicity (0=Swedish 1=2nd generation immigrants). 



Ethical approval 

The study outline was approved by the Ethics committee 
at the University of Lund in 1994 and 1998 and by The 
Regional Ethical Review Board in Linkoping 2007. 

Results 

The frequencies of 5-HTTLPR genotypes were 1/1 31%, 
s/1 47% and s/s 22%. The frequencies of BDNF Val66Met 
genotypes were Val/Val 65%, Val/Met 31%, Met/Met 4%. 

Bivariate analysis 

In bivariate analysis, maternal symptoms of depression 
and anxiety at the 12-year follow up increased the risk for 
internalizing problems in the children more than fivefold 
(OR 5.92, CI 3.59-9.77). There was an increase of risk also 
for the subscales anxious depressed and withdrawn de- 
pressed (OR 4.76, CI 2.95-7.70; OR 4.54, CI 2.82-7.30). A 
fivefold risk increase was seen for externalizing problems 
(OR 5.20, CI 3.26-8.28) (Table 1). 

Using the screening cut off of 10 on EPDS no effect was 
seen for maternal depressive symptoms postpartum on in- 
ternalizing problems and the subscale anxious depressed 
in the children. There was an increase of risk for with- 
drawn depressive symptoms (OR 2.04, CI 1.09-3.81) and 
externalizing symptoms (OR 2.07, CI 1.13-3.80) if the 
mother reported on depressive symptoms postpartum 
(Table 1). Using the higher cut off of 13 reduced the num- 
ber of EPDS positive mothers from 92 to 35. The associ- 
ation between maternal depressive symptoms postpartum 
and all CBCL scales were strengthened using the clinical 
cut off of 13 on the EPDS (Table 1). 

Analysis further showed that 5-HTTLPR s/s carriers 
had more internalizing problems at age 12 compared to 
1/1 carriers (OR 4.43, CI 2.09-9.38). The same applied for 
the subscales anxious depressed and withdrawn de- 
pressed (OR 2.32, CI 1.22-4.42; OR 2.48, CI 1.29-4.77) as 
well as for externalizing problems (OR 2.4, CI 1.23-4.75). 
Differences between s/1 and 1/1 carriers were not signifi- 
cant. No effect was seen for BDNF Val66Met genotype 
on child behavior problems (Table 1). 



Traumatic life events increased the risk for internalizing 
and externalizing problems in the children (OR 1.70, CI 
1.04-2.80; OR 1.62, CI 1.02-2.56). For the subscale anxious 
depressed, there was an increase of risk (OR 1.72, CI 1.07- 
2.72), but not concerning the subscale withdrawn de- 
pressed (Table 1). 

Maternal life stress at baseline increased the risk for in- 
ternalizing problems in the children (OR 4.03, CI 1.95- 
8.35). A risk increase was also seen for the subscale anxious 
depressed (OR 3.53, CI 1.71-7.28), but not for withdrawn 
depressed and externalizing problems (Table 1). 

Children with a non-Swedish background had an in- 
creased risk for both internalizing and externalizing prob- 
lems (OR 1.93, CI 1.02-3.67; OR 1.89, CI 1.04-3.46) 
compared to children with a Swedish background. They 
also indicated an increased risk on the subscale withdrawn 
depressed (OR 2.22, CI 1.22-4.01) (Table 1). Boys had an 
increased risk for externalizing problems compared to 
girls (OR 1.86, CI 1.17-2.95), whereas girls had an in- 
creased risk on the subscale anxious depressed (OR 0.59, 
CI 0.37-0.96) (Table 1). 

Multivariate analysis 

In the stepwise multivariate analysis using cut off of 10 on 
EPDS, maternal symptoms of depression and anxiety at 
the child's age of 12 increased the risk for both internaliz- 
ing and externalizing behavior problems in the children 
(OR 5.72, CI 3.30-9.91; OR 5.47, CI 3.40-8.78). When 
dividing the internalizing scale into the subscales anxious 
depressed and withdrawn depressed, the effect remained 
only for the latter (OR 4.27, CI 2.57-7.10) (Figure 1). 

5-HTTLPR s/s carriers had a more than 4-fold increased 
risk of internalizing problems at age 12 compared to 1/1 
carriers (OR 4.73, CI 2.14-10.48). Likewise, a risk increase 
was seen for the subscales anxious depressed and with- 
drawn depressed (OR 2.41, CI 1.25-4.62; OR 2.30, CI 1.18- 
4.51) (Figure 1). Differences between s/1 and 1/1 carriers 
were not significant. The BDNF Val66Met polymorphism 
did not show an impact on child behavior, and was not 
included in the final models. 

Maternal life stress at baseline increased the risk for 
internalizing problems in 12-year old children (OR 3.42, 
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CI 1.43-8.21). This was true also for the subscale anxious 
depressed (OR 3.91, CI 1.81-8.44) (Figure 1). No effect 
was seen for LITE on any of the models. 

Boys had an increased risk of externalizing problems 
compared to girls (OR 1.84, CI 1.14-2.97). No other dif- 
ferences relating to sex were detected. 

The only change when using a cut off of 13 on EPDS 
was seen concerning the outcome in the model for the 
subscale anxious depressed. The significance of 5- 
HTTLPR disappeared while both being girl (OR 0.53, CI 
0.32-0.87) and maternal symptoms of depression and anx- 
iety at the child's age of 12 (OR 4.62 CI 2.81-7.58) in- 
creased the risk for being anxious depressed according to 
the CBCL subscale. No separate significant effect was seen 
for EPDS using either of the two cut offs. 

There was no evidence for gene-by-environment interac- 
tions, or gene-by-gene-by-environment interactions. Fur- 
thermore, we tested whether there was an interaction 
between ethnic background and 5-HTTLPR and BDNF 
Val66Met polymorphisms respectively on the outcome of 
child behavior problems, but no such interaction was found 
(for all factors and interactions tested, see Additional file 1). 

Discussion 

This study used data from a longitudinal birth cohort 
study in order to examine predictors of internalizing and 
externalizing symptoms in 12-year old children. Environ- 
mental as well as genetic factors were included in order 
to test for gene-by-environment and gene-by-gene-by- 
environment interactions. The results of the study can 
be summarized in three main findings. 

First and foremost, maternal symptoms of depression 
and anxiety predicted child behavior problems at age 12. 
Numerous researchers have been able to show the import- 
ance of maternal mental health for child development and 
wellbeing. Our results indicate that in the present popula- 
tion, maternal symptoms of postpartum depression (ana- 
lyzed with different cut off scores on the EPDS) did not 
have a long-term impact on child behavior, i.e., no increase 
of risk was seen at age 12. Symptoms of depression in 
mothers at the 12-year follow-up could possibly reflect re- 
current or chronic depressive problems, which would put 
the child under greater stress than would be the case with 
a single episode depression. In clinical practice, we suggest 
that besides focusing on the child's mental health, it is 
important also to ask about depressive symptoms in 
the mother. 

Second, we found a main gene effect of 5-HTTLPR on 
child behavior at age 12. Since the finding of an inter- 
action between 5-HTTLPR and stressful life events [17], 
focus of research has shifted to gene-by-environment 
interaction effects. The present finding does have sup- 
port in the literature [16,31-34], however most gene-by 
-environment studies have not shown a main effect of 5- 



HTTLPR on depression [20,21]. Recently, Duncan and 
Keller reviewed the literature on candidate gene-by-en- 
vironment interaction (cGxE) studies in psychiatry over 
a decade (2000-2009), and based on the analyses of 103 
studies suggested that "well-powered direct replications 
deserve more attention than novel cGxE findings and 
indirect replications" [58]. The current study attempted 
testing a previous hypothesis and replicating previous 
findings in a large and longitudinal population-based in- 
dependent study/sample. The set-up of the current study 
can help to understand the generalizability of previous 
findings since it used a virtually similar set-up compared 
to previous studies with regard to phenotypic variables, 
genetic polymorphisms, statistical model, environmental 
moderator, and inclusion of both sexes [58]. This repre- 
sents a major strength of the current study which might 
contribute to the drawing of clearer conclusions in the 
context of future meta-analytical studies. The present 
study did not find any GxE or GxGxE effects but rather 
a main effect of 5-HTTLPR on internalizing symptoms 
as for example showed by Hoefgren et al. in a case-con- 
trol study on depression [31]. 

Third, traumatic life events showed an impact on child 
behavior in bivariate models, but not in the multivariate 
analyses, indicating that in this material, life events did not 
explain as much of the variance as other factors. We used 
the mothers' reports on traumatic life events experienced 
by the children. Previous studies have shown differences 
between parents' reports and children's self-reports on ex- 
posure to violence [59]. Considering this, self-reports 
might have contributed to a more fair estimation of trau- 
matic life events experienced by the child. Given the 
known relationship between life events and depressive 
symptoms, one would have expected an effect on the in- 
ternalizing scales [12,15,17]. 

We included measure of externalizing problems since 
previous research has shown high comorbidity between de- 
pressive and externalizing disorders in adolescents [60] and 
since boys more often express their psychological problems 
or stress by showing externalizing problems [2]. Further- 
more, Laucht et al [61] reported that adolescents with de- 
pressive disorders were more likely to have a history of 
externalizing problems. As expected, the analysis showed 
that boys expressed significantly more externalizing prob- 
lems than girls, but there were no differences relating to 
sex for internalizing problems. 

Despite the dropout level, the number of its participants 
and the longitudinal approach strengthens the study. There 
are some limitations that need to be mentioned. Although 
dropout rates of the same magnitude in longitudinal, multi 
wave studies are common [7], the skewed dropout in this 
study is a limitation - psychosocially disadvantaged families 
and families of foreign origin were less likely to participate 
in the follow-up. Considering the known connection 



Agnafors et al. Child and Adolescent Psychiatry and Mental Health 201 3, 7:1 0 
http://www.capmh.eom/content/7/1/10 



Page 7 of 9 



a. Internalizing 



HSCL-25t 
5.72 (95% CI=3.30-9.91) 



5-HTTLPR $ 
4.73 (95%CI=2.14-10.48) 



LSSt 

3.42 (95% CI=1.43-8.21) 




b. Anxious depressed 



5-HTTLPR $ 
2.41 (95% CI=1.25-4.62) 



LSSt 

3.91 (95% CI=1.81-8.44) 




Anxious 

depressed 
r 2 =0.05 



V 



C. Withdrawn depressed 



d. Externalizing 



5-HTTLPR $ 
2.30 (95% CI=1. 18-4.51) 



HSCL-25t 
4.27 (95% CI=2.57- 



7.10) 



HSCL-25^ Sex * 
5.47 (95% CI=3.40-8.78) [ 1.84(95% CI=1.14-2.97) 



Withdrawn 

depressed 





Figure 1 Odds ratios for the studied CBCL outcomes and corresponding measure of r2. $Comparing s/s with l/l as reference. tComparing 
risk with no risk as reference. "Comparing boys with girls as reference. Note: CBCL = Child Behavior Checklist, r 2 = the amount of variance explained 
by the model, 5-HTTLPR = serotonin transporter gene-linked polymorphic region, HSCL-25 = Hopkins Symptom Checklist 25, LSS = Life Stress Score. 
Multivariate analysis. CBCL scales presented are a) internalizing symptoms, internalizing subscales b) anxious depressed and c) withdrawn depressed, 
and d) externalizing symptoms. 



between psychosocial strain and behavior problems, it is 
plausible to think that a more representative participating 
population would have strengthened the results. 

In general, both child and parent-reports indicated low 
frequencies of behavior problems. Reports on good men- 
tal health in 12-year old children have been published 
previously. Costello et al. [2] examined mental health in 
9-16 year old children using parents reports on child 
behavior (CBCL) and clinical interviews, and found the 
lowest prevalence of depression in 12-year olds (0.4%). 
Conducting the follow-up at ages eight-ten or during 
adolescence, where previous reports have shown more 
behavior problems, might have increased the prevalence 
rate and thereby rendering a larger group for different 
analyses. 

Since there are studies indicating that depressed mothers 
are prone to overestimate behavior problems in their chil- 
dren [62], another limitation of the present study is that 
mothers' reports, not self-reports or teachers' reports, on 
child behavior and traumatic life events were used. We 



decided to use mothers as informants since they have 
known their children for 12 years around the clock com- 
pared to the current teachers that have known the children 
for about just over a year 5-6 hours a day, five days a week. 
Another support for this strategy was that although the as- 
sociation between maternal depressive symptoms and a 
positive mother-child reporting discrepancy has been 
demonstrated, in the general population this association is 
small and does not bias research, using maternal reports 
on child problems [62]. The chosen methodology used in 
the study also allowed following data from the same infor- 
mants from baseline to the 12-year follow-up. 

Feder, Nestler and Charneys report of 2009 called 
for an increased understanding of the psycho-biological 
factors which are involved in risk and resiliency for 
psychiatric disorders [63], and Duncan and Kellers re- 
view called for well-powered direct replications of cGxE 
findings [58]. Thus, many studies are still to be carried 
out to shed light on the psychogenetic underpinning of 
depression. 
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Conclusion 

In conclusion the present study suggests that concurrent 
maternal symptoms of depression and anxiety are an im- 
portant risk for behavior problems in children, which needs 
to be taken into account in clinical practice. Furthermore, 
we found a main effect of 5-HTTLPR on internalizing 
symptoms in 12 year old children, a finding that need to be 
confirmed in future studies. 
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Additional file 1: Factors and interactions included in multivariate 
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between factors tested in multivariate analysis. 
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